Phase transitions in the mesoscopic superconducting square 
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We solve the Ginzburg-Landau equation (GLE) for the mesoscopic thin film of the square shape 
in the magnetic field. In the limit of Ginzburg-Landau parameter k — > oo we find a series of first 
and second order phase transitions as temperature and/or magnetic field changes. First order phase 
transitions between giant flux states can be described with a simple variational procedure. We 
discuss the similarity with rotating liquid He 4 and derive a simple formula for H c \. We identify 
order parameters based on symmetry arguments and we propose Landau functional describing the 
second order phase transition. 



PACS: 74.60. Ec 74.25. Ha 74.80.-g 

Advances in nanotcchnology and constantly shrink- 
ing semiconductor devices have motivated researches to 
study properties of mesoscopic superconducting samples. 
One line of research in this field has focused on the prob- 
lem of the phase transitions in the mesoscopic supercon- 
ducting sample under the influence of the external mag- 
netic field [0 . There are two characteristic limits in which 
phase transitions have different properties. If the size of 
the sample a ^> £, where £ is the superconducting coher- 
ence length, and applied field are large enough, there are 
many vortices in the sample. In this case long-range in- 
teraction between vortices and image vortices is screened 
by spontaneous creation of vortex loops near the sample 
boundary 0] . It leads to the decrease of the surface bar- 
rier for the vortex to penetrate into the sample. In the 
opposite case, when a ~ £ <C A, with A being the Lon- 
don penetration depth, there are only few vortices in the 
sample. The standard Abrikosov approach || must be 
modified because of the strong influence of the sample 
boundaries. In this case magnetic field is almost uni- 
form throughout the sample and long-range forces be- 
tween vortices are not important. Thermodynamics of 
this system is determined by the short-range repulsion of 
vortices and Bean-Livingston barrier forces B . 

Different approaches have been applied for the inves- 
tigation of phase transitions in the latter limit. Most 
of them consider disk geometry. Buzdin and Brison ap- 
plied electrostatic formalism to consider influence of the 
barrier on the vortex structure of the thin superconduct- 
ing disk J5|. Within this approach vortices are replaced 
by the hard-core particles interacting through Coulomb 
forces and the giant vortex has never been discussed. Nu- 
merical solution of GLE for the same geometry reveals a 
series of the first and second order phase transitions in 
the superconducting disk. Such transitions take place 
between giant vortex states with different vorticity as 
well as between a giant vortex state and a multi-vortex 
state as the external field changes We emphasize 

an important difference between the disk and the square 
geometry. Solution of the linearized GLE, describing the 



nucleation of superconducting order parameter near H C 2 
line for the disk, always corresponds to the giant vortex 
state. On the other hand, as it was demonstrated by 
Chibotaru et al. jjj , there are many well separated zeros 
of the order parameter in the case of the square sample. 
Consequently, the behavior of the square sample near 
H C 2 line should be qualitatively different from the disks. 
On the basis of the solution of the linearized GLE the 
appearance of the antivortex in the center of the sample 
has been predicted 0. 

In the this paper we investigate phase transitions in 
the superconducting film of a square shape as a function 
of temperature T and external magnetic field H in the 
limit a ~ ( C A. We solve GLE for the thin super- 
conducting film with the thickness d -C £. We also as- 
sume that superconducting film is an extreme type two 
superconductor, where the Ginzburg-Landau parameter 
k = A/£ 3> 1. We show that configuration with one anti- 
vortex in the center and four vortices on the diagonals of 
the square is unstable when we move away from the H C 2 
line and nonlinear term in the GLE is considered. On the 
contrary, at higher magnetic field, the configuration with 
four vortices on diagonals of the square remains stable. 
We find a sequence of phase transitions of the first order 
between giant vortex states as well as between multi- 
vortex states with different vorticity. Second order phase 
transition takes place when a giant vortex state splits into 
a multi- vortex state with simultaneously breaking the C4 
symmetry. Such transitions are discussed in terms of the 
phcnomenological theory of Landau. 

GLE for dimensionless complex order parameter ip has 
the following form: 

9tt A 

here £ = 4t „ | a | , ol is temperature dependent parameter 
of Ginzburg-Landau expansion for the free energy, $0 is 
the flux quantum, A is the vector potential H = V x A. 
The second GLE equation for the vector potential can be 
written as: 
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Since we consider the case of a small mesoscopic square 
where a ~ £ <C A, the magnetic field is uniform in the 
film. The correction to the external field is of the order of 
1/k 2 and may be found by solving Eq.(l) while assuming 
uniform magnetic field and substituting the solution of 
Eq.(l) to Eq.(2). Such a solution is equivalent to the 
expansion of the free energy in 1/k 2 series. In addition 
to Eq.(l) we have to supply the boundary condition for 
the superconductor-insulator junction: 
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where n is normal vector to the surface of the sample. 

Introducing N x N discrete points in the square 
we rewrite Eq.(l) in the form of nonlinear discrete 
Schrodinger equation: 
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where summation index 1 = (±1,0), (0, ±1) points to- 
ward nearest neighbours and 
ti lt i = (£N/a) 2 exp(-^ J} 1 A(r)dr) |]. Equivalent dis- 
cretization of the boundary conditions, Eq.(3), provides 
an additional equation which can be directly solved and 
substituted into Eq.(4). As a result, equation close to 
the boundary are slightly different from the 'bulk': 
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- e(i)ti iV>i - ipi + ^ilV'il 2 = 0, 
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where ipi = if i is outside of the sample, e(i) = 4— 6^,1 — 
Si a ,N—8i y ,i—6i yi N andi = (i x = 1, . . . , N, i y = 1, . . . ,N). 
There is one important advantage of such a treatment of 
the boundary condition. When neglecting the nonlinear 
term in the Eq. (4) , the system of linear equations reduces 
to the problem of eigenvalues and eigenfunction of the 
hermitean matrix. On the other hand, the solution of 
nonlinear equations requires iterations and inversion of 
the hermitean matrix. 

Let us first discuss the solution of the linearized GLE 
and compare our results with previous studies J5J. The 
lowest eigenvalue of the linear GLE determines the upper 
critical field of the sample. We have calculated eigenval- 
ues of the linear problem expressed in units (a/£(T)) 2 as 
a function of the dimcnsionlcss external magnetic field 
h = <i>/$o where $ is the total flux through the sample. 
Our results for a few lowest eigenvalues agree within the 
linewidth with the results of Ref . jl| . Spatial pattern of 
the order parameter is also similar. For the field h ~ 5.5 
we have observed 5 zeros of the order parameter near the 
center of the square. The solution corresponds to 4 vor- 
tices on the diagonals and one anti-vortex in the center 
of the square with total vorticity m = 3 H . The distance 
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FIG. 1. a) Difference in the free energy AF in units of 

<J> 2 a 2 d 

(47r)° a A^ 2 between solutions with different integral vorticities. 
For the case of h = 5.5, AF = F(m = 3) - F(m = 2) and for 
h = 7, AF = F(m = 4) — F(m = 3). b) Order parameter r\ vs. 
(°/0 2 f° r the solution with m = 2 calculated at h = 5.5 and 
two dimensional order parameter r] x and r\ y for the solution 
with m = 3 calculated at h = 7. c) The second derivative 
of the free energy A = d 2 F/d((a/^) 2 ) 2 for the solution with 
m = 2 at h = 5.5. 



between vortices is of the order of 5 ~ 0.12£ << £. All 
these zeros are pinned in the region of the square where 
\ip(x, y) | is small (four orders less then the value of the or- 
der parameter near the sample boundary). This indicates 
that vortex-anti- vortex structure becomes unstable when 
we move away from the H C 2 line and nonlinear term and 
1/k corrections are considered. Moreover, the value of 
the order parameter and the screening current are small 
near zeros. In that case correction to the external field is 
determined by the current flowing around all 5 zeros of 
the order parameter rather than by the current between 
them. Therefor, we do not expect suppression of the field 
in the core of the anti- vortex. 

In the following we consider the changes in the vortex 
structure when we move away from the H C 2 line. Tak- 
ing into account the nonlinear term, we were unable to 
detect more than one zero of the order parameter near 
the origin of the square for the value of the field h = 5.5. 
Consequently, the solution with more than one zero sur- 
vives only very close to H c2 line. We do not expect any 
phase transition at the point where all zeros are joined 
together since total flux and the symmetry of the solu- 
tion do not change. This situation is different from the 
case of higher field h = 7, where solution with four zeros 
of the order parameter survives far from the H C 2 line. In 
fact, the giant flux solution was not detected in that case. 

In the Fig. (Qa) we plot the difference in the free en- 
ergy between solutions with different integral vorticity 
AF = F(m = 3) - F(m = 2) as a function of (a/f) 2 for 
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the fixed magnetic field h = 5.5. As it is clearly seen from 
the figure, near (a/£) 2 ~ 43 a first order phase transition 
takes place. At that point the high-temperature phase 
corresponding to the giant vortex with m = 3 becomes 
metastable, and the phase corresponding to giant vor- 
tex, shown in Fig. (|]top), with m — 2 becomes a ground 
state. At this point, the slope of the first derivative of 
the free energy as a function of T is discontinuous, which 
corresponds to the latent heat of the transition. With 
further decrease of the temperature and £, the second 
transition takes place. At (a/£) 2 ~ 66 the giant vor- 
tex located in the center splits along one of the diago- 
nals, Fig. (H bottom). This transition is the second order 
phase transition, where \ip\ 2 is no longer invariant un- 
der four-fold axis of the square. The phase transition 
is clearly observed by computing the order parameter 
7] = J xy\ip(x,y)\ 2 dxdy, presented in Fig. (0b). A non- 
vanishing rj is followed by a jump in the second deriva- 
tive of the free energy A = d 2 F '/ d((a / £) 2 ) 2 , presented in 
Fig. (|l|c). We estimated the magnitude of the effect in 

terms of the specific heat jump: AC/a 2 d = , 

where T c is the critical temperature of the transition de- 
termined by the condition (a/£) 2 ~ 66, and T c is the 
critical temperature of the sample at h = 0. Landau 
functional, describing this phase transition, is defined as 
F = a\T] 2 -\- (5\r\ , where a.\ and f3\ are Landau coefficients 
with ot\ oc (T — T c ). This phase transition corresponds 
to the one dimensional corepresentation B of the nonuni- 
tary C 4t ,(C 4 ) group. 

Since we restrict our calculation to the limit of n —> oo 
and A » a, long range forces (r 3> £) between vortices 
are irrelevant. At short distances r ~ £ there is short 
range repulsion between them. There is another interac- 
tion of vortices with the boundaries of the sample, known 
as Bean-Livingston barrier Q . In the vicinity of H C 2 line 
the interaction with the boundaries is larger than the 
repulsion between vortices. This is the case for the mag- 
netic field h — 5.5, where the giant vortex with m = 3 
is located in the center of the sample. With the decrease 
of £ the energy difference between two different solutions 
with different vorticity decreases and the first order phase 
transition to a giant flux state with vorticity m = 2 takes 
place at (a/£) 2 — 43. A giant flux state (with lower vor- 
ticity) remains stable, since the interaction with the sam- 
ple boundaries still dominates over the repulsion between 
vortices. Further decrease of £ leads to the decrease of 
the interaction of the vortices with the boundaries. Bean- 
Livingston force decreases and repulsion of the vortices 
in the giant flux state starts to dominate. As a result of 
the interplay of the vortex-vortex repulsion and the re- 
pulsion of vortex from the boundaries, the second order 
phase transition (at (a/£) 2 — 66) between the giant vor- 
tex state and the multivortex state takes place, preserv- 
ing the integral vorticity. Separation between vortices is 
determined by the vortex-vortex repulsion that tends to 




FIG. 2. \ip(x,y)\ calculated at fixed magnetic field h = 5.5. 
Presented are solutions with lowest free energy calculated at 
(a/£) 2 = 50 with one giant flux m = 2 (top) and (a/f;) 2 = 90 
with two separated fluxes each carrying m = 1 (bottom). 
Contours represent \tj)(x,y)\ = 0.1, . . . ,0.9. 

separate them as far as possible, while in the contrary, 
repulsion from the boundaries prevents vortices from ap- 
proaching the boundaries. 

Situation is different when external magnetic field is in- 
creased to h ~ 7. In that case near H c2 line ground state 
corresponds to a multivortex state with the total vortic- 
ity m = 4, Fig. (||top). When temperature decreases, the 
first order phase transition takes place at (a/£) 2 = HO, 
see Fig. (0a). At this point multivortex state with 
Til — 4 becomes unstable while the multivortex state with 
rn = 3, presented in Fig. (^ bottom), represents the solu- 
tion with the lowest free energy. Apart from the change 
of the vorticity, the symmetry is also reduced at the tran- 
sition point. Phase transition takes place in accordance 
with joint corepresentation of the nonunitary group 
Ci v {Ci). Consequently, four orientations of the pseu- 
dodipolar moment of the vortices are possible. Two com- 
ponent order parameter corresponding to given change of 
symmetry, presented in Fig. (^b), can be determined as 
follows: i]x = J x\ijj(x,y)\ 2 dxdy, r) y = / y\ip(x, y)\ 2 dxdy. 
Free energy in that case depends on the vorticity m and 
order parameter r\ Xl r\ y . F(m — A,r] x ,r]y) always has a 
minimum at rj x — r\ y — 0. For m — 3, F(m — 3, r] x , rj y ) — 
ot2(vl + Vy) + foivt + Tl$) + l2vlvl, where a 2 , (h, li are 
Landau coefficients, a 2 oc (T — T c ) where T c is the tem- 
perature of the transition between the giant vortex state 
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FIG. 3. \if>(x, y)\ calculated at fixed magnetic field h = 7. 
Presented are solutions with lowest free energy calculated at 
( a /0 2 = 100 ( to P) and ( a /£) 2 = 120 (bottom). Contours are 
defined as in Fig. (Q). 

and the multivortex state for the case of m = 3. T c is 
determined by the condition (a/£) 2 ~ 56. This transition 
is unobservable because T c is lower than the transition 
temperature for the first order phase transition where 
the vorticity m changes form m — 3 to m = 2 (compare 
Figs. (0a) and (0b)). 

The first order phase transition between giant vortex 
states with different vorticity can be described qualita- 
tively on the basis of a simple variational function for 
the order parameter. Spatial dependence of the order 
parameter in the giant vortex state with vorticity m can 
be approximated by the function: 

Mr 6) = { (r/O m e*P (im<f>); r < f, 
y ' \ exp (zm</>) ; r > £. w 

Substituting this function into the Ginzburg-Landau 
functional and keeping only leading terms in (a/£) 2 , we 
obtain a simple expression for the free energy which 
should be minimized as a function of vorticity m: 

F ~ (m 2 In (a 2 /< 2 )/2 - $ . ro /$ ). (7) 

Minimization of the Eq.(7) in the large m 3> 1 limit 
provides the expression for the vorticity: 

m^/ln(a 2 /< 2 ). (8) 



Minimization of the Eq.(7) for h = 5.5 provides that 
the phase transition from m = 3 to m — 2 takes place at 
(°/£) 2 = 28 which is lower than calculated value (a/£) 2 — 
43. Nevertheless, this simple variational formula provides 
simple explanation of the transition with the change of 
vorticity in the giant vortex. It is interesting to note 
that due to logarithmic term, Eq.(7) for the free energy 
is similar to the free energy of the rotating superfluid 
liquid Q . This similarity appears because A > a and all 
integrals are cut at a, rather than at A. Eq.(7) yields also 
the estimate for H c i, the field at which the first vortex 
appears in the sample. Substituting m — 1 to the Eq.(l) 
and solving equation F(m = 1) = we obtain: 

#ci = |j-ln(a 2 M 2 ). (9) 

This expression is similar to the bulk H c i where A 2 is 
substituted to a 2 /ir. 

In conclusion we have solved GLE in the limit of k — > 
co for the thin square film. We have predicted a series 
of the first and second order phase transitions with the 
change of temperature and magnetic field. On the basis 
of the symmetry we constructed Landau functional for 
the second order transitions. First order transitions are 
described on the basis of variational estimates. 
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